In-class Exercise 4: Preparing Spatial Interaction Modelling Variables
1. Objective
- Perform geocoding data downloaded from data.gov.sg using SLA OneMap API,
- Convert an aspatial data into a simple feature tibble data.frame,
- Perform point-in-polygon count analysis,
- Append the propulsiveness and attractiveness variables onto a flow data, and
- Calibrate Geographically Weighted Poisson Regression
2. Getting Starting
httr used to communicate with web server.
3. Geocoding using SLA API
Next, the code chunk below will be used to combine both found and not_found data.frames into a single tibble data.frame called merged. At the same time, we will write merged and not_found tibble data.frames into csv file format for subseuent use.
4. Converting an aspatial data into a single feature tibble data.frame
4.1 Importing and tidying schools data
4.2 Converting an aspatial data into sf tibble data.frame
Next, you will convert schools tibble data.frame data into a simple feature tibble data.frame called schools_sf by using values in latitude and longitude fields.
3 variables in output sf tibble data frame. Lon&lat combined into 1 point column.
4.3 Plotting a point simple feature layer
To ensure that schools sf tibble data.fame has been projected and converted correctly, plot the schools point data for visual inspection.
Reading layer `MPSZ-2019' from data source
`C:\kytjy\ISSS624\In-class_Ex\In-class_Ex4\data\geospatial'
using driver `ESRI Shapefile'
Simple feature collection with 332 features and 6 fields
Geometry type: MULTIPOLYGON
Dimension: XY
Bounding box: xmin: 103.6057 ymin: 1.158699 xmax: 104.0885 ymax: 1.470775
Geodetic CRS: WGS 844.4 Performing Point-in-Polygon Count Process
It is always a good practice to examine the summary statistics of the derived variable.
Min. 1st Qu. Median Mean 3rd Qu. Max.
0.000 0.000 0.000 1.054 2.000 12.000 The summary statistics above reveals that there are excessive 0 values in SCHOOL_COUNT field. If log() is going to use to transform this field, additional step is required to ensure that all 0 will be replaced with a value between 0 and 1 but not 0 neither 1.
5. Data Inetegration and Final Touch-Up
Reading layer `Business' from data source
`C:\kytjy\ISSS624\In-class_Ex\In-class_Ex4\data\geospatial'
using driver `ESRI Shapefile'
Simple feature collection with 6550 features and 3 fields
Geometry type: POINT
Dimension: XY
Bounding box: xmin: 3669.148 ymin: 25408.41 xmax: 47034.83 ymax: 50148.54
Projected CRS: SVY21 / Singapore TM6. Data Integration & Wrangling
7. Model Calibration
Rows: 14,734
Columns: 13
$ ORIGIN_SZ <chr> "AMSZ01", "AMSZ01", "AMSZ01", "AMSZ01", "AMSZ01", "AMS…
$ DESTIN_SZ <chr> "AMSZ01", "AMSZ02", "AMSZ03", "AMSZ04", "AMSZ05", "AMS…
$ MORNING_PEAK <dbl> 1998, 8289, 8971, 2252, 6136, 2148, 1620, 1925, 1773, …
$ dist <dbl> 50.0000, 810.4491, 1360.9294, 840.4432, 1076.7916, 805…
$ ORIGIN_AGE7_12 <dbl> 310, 310, 310, 310, 310, 310, 310, 310, 310, 310, 310,…
$ ORIGIN_AGE13_24 <dbl> 710, 710, 710, 710, 710, 710, 710, 710, 710, 710, 710,…
$ ORIGIN_AGE25_64 <dbl> 2780, 2780, 2780, 2780, 2780, 2780, 2780, 2780, 2780, …
$ DESTIN_AGE7_12 <dbl> 310.00, 1140.00, 1010.00, 980.00, 810.00, 1050.00, 420…
$ DESTIN_AGE13_24 <dbl> 710.00, 2770.00, 2650.00, 2000.00, 1920.00, 2390.00, 1…
$ DESTIN_AGE25_64 <dbl> 2780.00, 15700.00, 14240.00, 11320.00, 9650.00, 12460.…
$ SCHOOL_COUNT <dbl> 0.99, 2.00, 2.00, 1.00, 3.00, 2.00, 0.99, 0.99, 3.00, …
$ RETAIL_COUNT <dbl> 1.00, 0.99, 6.00, 0.99, 0.99, 0.99, 1.00, 117.00, 0.99…
$ geometry <LINESTRING [m]> LINESTRING (29501.77 39419...., LINESTRING …flow_data$FlowNoIntra <- ifelse(
flow_data$ORIGIN_SZ == flow_data$DESTIN_SZ,
0, flow_data$MORNING_PEAK)
flow_data$offset <- ifelse(
flow_data$ORIGIN_SZ == flow_data$DESTIN_SZ,
0.000001, 1)
inter_zonal_flow <- flow_data %>%
filter(FlowNoIntra >0)
inter_zonal_flow <- inter_zonal_flow %>%
rename(TRIPS = MORNING_PEAK,
DIST = dist)8. Origin (Production) Constrained SIM
orcSIM_Poisson <- glm(formula = TRIPS ~
ORIGIN_SZ +
log(SCHOOL_COUNT) +
log(RETAIL_COUNT) +
log(DIST) -1, # -1 to remove intersect for origin-constrained and destination-constrained
family = poisson(link="log"),
data= inter_zonal_flow,
na.action= na.exclude) # excludes any NAs in the data
summary(orcSIM_Poisson)
Call:
glm(formula = TRIPS ~ ORIGIN_SZ + log(SCHOOL_COUNT) + log(RETAIL_COUNT) +
log(DIST) - 1, family = poisson(link = "log"), data = inter_zonal_flow,
na.action = na.exclude)
Coefficients:
Estimate Std. Error z value Pr(>|z|)
ORIGIN_SZAMSZ01 19.8739840 0.0047627 4172.84 <2e-16 ***
ORIGIN_SZAMSZ02 20.5902203 0.0042786 4812.33 <2e-16 ***
ORIGIN_SZAMSZ03 20.2327026 0.0045531 4443.70 <2e-16 ***
ORIGIN_SZAMSZ04 19.7744438 0.0049837 3967.79 <2e-16 ***
ORIGIN_SZAMSZ05 19.6574529 0.0056396 3485.61 <2e-16 ***
ORIGIN_SZAMSZ06 19.9659115 0.0048946 4079.16 <2e-16 ***
ORIGIN_SZAMSZ07 18.6746164 0.0096316 1938.90 <2e-16 ***
ORIGIN_SZAMSZ08 19.2701601 0.0090776 2122.82 <2e-16 ***
ORIGIN_SZAMSZ09 19.9889467 0.0052858 3781.64 <2e-16 ***
ORIGIN_SZAMSZ10 20.3422035 0.0045778 4443.62 <2e-16 ***
ORIGIN_SZAMSZ11 18.3944113 0.0129212 1423.58 <2e-16 ***
ORIGIN_SZAMSZ12 18.3484209 0.0109652 1673.33 <2e-16 ***
ORIGIN_SZBDSZ01 20.9668587 0.0043388 4832.36 <2e-16 ***
ORIGIN_SZBDSZ02 20.4059518 0.0050601 4032.75 <2e-16 ***
ORIGIN_SZBDSZ03 20.6725514 0.0045276 4565.93 <2e-16 ***
ORIGIN_SZBDSZ04 21.6703853 0.0038930 5566.44 <2e-16 ***
ORIGIN_SZBDSZ05 20.7497445 0.0046085 4502.46 <2e-16 ***
ORIGIN_SZBDSZ06 20.9119361 0.0046432 4503.77 <2e-16 ***
ORIGIN_SZBDSZ07 18.9749815 0.0097896 1938.28 <2e-16 ***
ORIGIN_SZBDSZ08 19.1933901 0.0091312 2101.95 <2e-16 ***
ORIGIN_SZBKSZ01 19.5422606 0.0064732 3018.96 <2e-16 ***
ORIGIN_SZBKSZ02 20.1748913 0.0050076 4028.89 <2e-16 ***
ORIGIN_SZBKSZ03 20.3984624 0.0047226 4319.35 <2e-16 ***
ORIGIN_SZBKSZ04 19.6182212 0.0059652 3288.76 <2e-16 ***
ORIGIN_SZBKSZ05 19.6033818 0.0063181 3102.74 <2e-16 ***
ORIGIN_SZBKSZ06 19.7145224 0.0056372 3497.20 <2e-16 ***
ORIGIN_SZBKSZ07 20.4237448 0.0041912 4873.03 <2e-16 ***
ORIGIN_SZBKSZ08 19.7992538 0.0050405 3928.02 <2e-16 ***
ORIGIN_SZBKSZ09 19.7821586 0.0055558 3560.66 <2e-16 ***
ORIGIN_SZBLSZ01 17.7977276 0.0149058 1194.01 <2e-16 ***
ORIGIN_SZBLSZ02 17.4287491 0.0192364 906.03 <2e-16 ***
ORIGIN_SZBLSZ03 16.5884288 0.0459848 360.74 <2e-16 ***
ORIGIN_SZBLSZ04 17.7851626 0.0232823 763.89 <2e-16 ***
ORIGIN_SZBMSZ01 20.0751840 0.0052887 3795.89 <2e-16 ***
ORIGIN_SZBMSZ02 18.6956140 0.0066656 2804.80 <2e-16 ***
ORIGIN_SZBMSZ03 19.3204425 0.0054755 3528.56 <2e-16 ***
ORIGIN_SZBMSZ04 19.4724220 0.0049390 3942.59 <2e-16 ***
ORIGIN_SZBMSZ05 16.9581801 0.0168804 1004.61 <2e-16 ***
ORIGIN_SZBMSZ06 16.9898638 0.0181852 934.27 <2e-16 ***
ORIGIN_SZBMSZ07 19.2868403 0.0056231 3429.91 <2e-16 ***
ORIGIN_SZBMSZ08 19.1477543 0.0055918 3424.28 <2e-16 ***
ORIGIN_SZBMSZ09 18.7564539 0.0086298 2173.46 <2e-16 ***
ORIGIN_SZBMSZ10 18.3617854 0.0089250 2057.35 <2e-16 ***
ORIGIN_SZBMSZ11 18.9167941 0.0063340 2986.54 <2e-16 ***
ORIGIN_SZBMSZ12 18.7874661 0.0093024 2019.63 <2e-16 ***
ORIGIN_SZBMSZ13 19.5654046 0.0057517 3401.70 <2e-16 ***
ORIGIN_SZBMSZ14 19.0685619 0.0063346 3010.24 <2e-16 ***
ORIGIN_SZBMSZ15 19.4403124 0.0058147 3343.30 <2e-16 ***
ORIGIN_SZBMSZ16 18.4469203 0.0092638 1991.28 <2e-16 ***
ORIGIN_SZBMSZ17 18.3430175 0.0157692 1163.22 <2e-16 ***
ORIGIN_SZBPSZ01 20.1806714 0.0053660 3760.81 <2e-16 ***
ORIGIN_SZBPSZ02 19.8116707 0.0061485 3222.19 <2e-16 ***
ORIGIN_SZBPSZ03 19.8467602 0.0059769 3320.57 <2e-16 ***
ORIGIN_SZBPSZ04 20.4613200 0.0048398 4227.72 <2e-16 ***
ORIGIN_SZBPSZ05 20.5379711 0.0043769 4692.39 <2e-16 ***
ORIGIN_SZBPSZ06 18.8948034 0.0093668 2017.21 <2e-16 ***
ORIGIN_SZBPSZ07 19.4104568 0.0087961 2206.70 <2e-16 ***
ORIGIN_SZBSSZ01 20.0139503 0.0056561 3538.45 <2e-16 ***
ORIGIN_SZBSSZ02 20.2543885 0.0047198 4291.38 <2e-16 ***
ORIGIN_SZBSSZ03 19.5428803 0.0052713 3707.41 <2e-16 ***
ORIGIN_SZBTSZ01 20.0198045 0.0058541 3419.77 <2e-16 ***
ORIGIN_SZBTSZ02 19.3618525 0.0081472 2376.51 <2e-16 ***
ORIGIN_SZBTSZ03 19.5883853 0.0068935 2841.59 <2e-16 ***
ORIGIN_SZBTSZ04 18.7720238 0.0103909 1806.58 <2e-16 ***
ORIGIN_SZBTSZ05 18.8069026 0.0120628 1559.08 <2e-16 ***
ORIGIN_SZBTSZ06 18.7068633 0.0094575 1978.00 <2e-16 ***
ORIGIN_SZBTSZ07 17.6292257 0.0141551 1245.43 <2e-16 ***
ORIGIN_SZBTSZ08 18.6989374 0.0109610 1705.94 <2e-16 ***
ORIGIN_SZCBSZ01 18.2189868 0.0548317 332.27 <2e-16 ***
ORIGIN_SZCCSZ01 18.9734563 0.0139450 1360.59 <2e-16 ***
ORIGIN_SZCHSZ01 19.5955119 0.0121035 1619.00 <2e-16 ***
ORIGIN_SZCHSZ02 19.3320960 0.0081620 2368.55 <2e-16 ***
ORIGIN_SZCHSZ03 21.2164518 0.0063552 3338.43 <2e-16 ***
ORIGIN_SZCKSZ01 20.1046845 0.0049333 4075.29 <2e-16 ***
ORIGIN_SZCKSZ02 20.5371946 0.0050256 4086.53 <2e-16 ***
ORIGIN_SZCKSZ03 20.7210560 0.0042184 4912.07 <2e-16 ***
ORIGIN_SZCKSZ04 21.4013886 0.0042524 5032.80 <2e-16 ***
ORIGIN_SZCKSZ05 20.9413146 0.0049434 4236.18 <2e-16 ***
ORIGIN_SZCKSZ06 20.2557727 0.0071832 2819.88 <2e-16 ***
ORIGIN_SZCLSZ01 19.3383703 0.0076634 2523.46 <2e-16 ***
ORIGIN_SZCLSZ02 18.5226956 0.0135522 1366.77 <2e-16 ***
ORIGIN_SZCLSZ03 19.0225512 0.0080145 2373.51 <2e-16 ***
ORIGIN_SZCLSZ04 20.7981505 0.0042400 4905.22 <2e-16 ***
ORIGIN_SZCLSZ05 18.3015625 0.0146815 1246.58 <2e-16 ***
ORIGIN_SZCLSZ06 20.8207386 0.0039567 5262.09 <2e-16 ***
ORIGIN_SZCLSZ07 19.6728958 0.0054199 3629.76 <2e-16 ***
ORIGIN_SZCLSZ08 20.0851929 0.0056956 3526.43 <2e-16 ***
ORIGIN_SZCLSZ09 18.5749589 0.0165415 1122.93 <2e-16 ***
ORIGIN_SZDTSZ02 15.8276209 0.0833992 189.78 <2e-16 ***
ORIGIN_SZDTSZ03 16.2512838 0.0737972 220.22 <2e-16 ***
ORIGIN_SZDTSZ13 16.7744385 0.0312450 536.87 <2e-16 ***
ORIGIN_SZGLSZ01 18.2368248 0.0096104 1897.62 <2e-16 ***
ORIGIN_SZGLSZ02 19.8705255 0.0049014 4054.06 <2e-16 ***
ORIGIN_SZGLSZ03 19.8249435 0.0053109 3732.85 <2e-16 ***
ORIGIN_SZGLSZ04 20.7800335 0.0041261 5036.20 <2e-16 ***
ORIGIN_SZGLSZ05 20.6040494 0.0043049 4786.23 <2e-16 ***
ORIGIN_SZHGSZ01 20.0273475 0.0044824 4468.04 <2e-16 ***
ORIGIN_SZHGSZ02 20.2480656 0.0044575 4542.47 <2e-16 ***
ORIGIN_SZHGSZ03 20.0756442 0.0049003 4096.81 <2e-16 ***
ORIGIN_SZHGSZ04 20.7577748 0.0040465 5129.84 <2e-16 ***
ORIGIN_SZHGSZ05 20.9779992 0.0040123 5228.42 <2e-16 ***
ORIGIN_SZHGSZ06 19.7403058 0.0054229 3640.20 <2e-16 ***
ORIGIN_SZHGSZ07 20.1896268 0.0046051 4384.22 <2e-16 ***
ORIGIN_SZHGSZ08 19.8646492 0.0052403 3790.72 <2e-16 ***
ORIGIN_SZHGSZ09 18.3647736 0.0069196 2654.04 <2e-16 ***
ORIGIN_SZHGSZ10 16.8720475 0.0421046 400.72 <2e-16 ***
ORIGIN_SZJESZ01 20.2673794 0.0046723 4337.79 <2e-16 ***
ORIGIN_SZJESZ02 20.0595982 0.0046503 4313.61 <2e-16 ***
ORIGIN_SZJESZ03 19.9128778 0.0049848 3994.75 <2e-16 ***
ORIGIN_SZJESZ04 18.5053667 0.0099227 1864.94 <2e-16 ***
ORIGIN_SZJESZ05 17.8172930 0.0138840 1283.29 <2e-16 ***
ORIGIN_SZJESZ06 20.0124157 0.0045009 4446.36 <2e-16 ***
ORIGIN_SZJESZ07 18.1821423 0.0117267 1550.49 <2e-16 ***
ORIGIN_SZJESZ08 18.8713046 0.0116456 1620.46 <2e-16 ***
ORIGIN_SZJESZ09 20.5535527 0.0048456 4241.72 <2e-16 ***
ORIGIN_SZJESZ10 18.4922322 0.0191243 966.95 <2e-16 ***
ORIGIN_SZJESZ11 18.2891211 0.0197114 927.85 <2e-16 ***
ORIGIN_SZJWSZ01 20.4912737 0.0063102 3247.35 <2e-16 ***
ORIGIN_SZJWSZ02 20.8236694 0.0042249 4928.82 <2e-16 ***
ORIGIN_SZJWSZ03 21.2587613 0.0039733 5350.40 <2e-16 ***
ORIGIN_SZJWSZ04 20.3816464 0.0046199 4411.67 <2e-16 ***
ORIGIN_SZJWSZ05 18.0607448 0.0128857 1401.61 <2e-16 ***
ORIGIN_SZJWSZ06 18.7015202 0.0107614 1737.83 <2e-16 ***
ORIGIN_SZJWSZ07 17.3991822 0.0277096 627.91 <2e-16 ***
ORIGIN_SZJWSZ08 21.8044465 0.0037356 5836.95 <2e-16 ***
ORIGIN_SZJWSZ09 21.5414930 0.0036033 5978.19 <2e-16 ***
ORIGIN_SZKLSZ01 20.0307712 0.0047868 4184.59 <2e-16 ***
ORIGIN_SZKLSZ02 19.0634769 0.0062318 3059.05 <2e-16 ***
ORIGIN_SZKLSZ03 19.2685700 0.0057172 3370.25 <2e-16 ***
ORIGIN_SZKLSZ04 17.7085067 0.0119809 1478.06 <2e-16 ***
ORIGIN_SZKLSZ05 18.6384471 0.0107596 1732.26 <2e-16 ***
ORIGIN_SZKLSZ06 13.7280296 0.1857160 73.92 <2e-16 ***
ORIGIN_SZKLSZ07 18.6425146 0.0084952 2194.47 <2e-16 ***
ORIGIN_SZKLSZ08 18.0928506 0.0101567 1781.37 <2e-16 ***
ORIGIN_SZLKSZ01 17.8907138 0.0397083 450.55 <2e-16 ***
ORIGIN_SZMDSZ01 18.7605188 0.0285455 657.22 <2e-16 ***
ORIGIN_SZMDSZ02 19.1533927 0.0102815 1862.90 <2e-16 ***
ORIGIN_SZMDSZ03 17.8404982 0.0169690 1051.36 <2e-16 ***
ORIGIN_SZMPSZ01 19.0765941 0.0083937 2272.74 <2e-16 ***
ORIGIN_SZMPSZ02 19.2162527 0.0068331 2812.24 <2e-16 ***
ORIGIN_SZMPSZ03 19.9965344 0.0054569 3664.44 <2e-16 ***
ORIGIN_SZMUSZ02 15.9130765 0.1037472 153.38 <2e-16 ***
ORIGIN_SZNTSZ01 17.0840999 0.0352513 484.64 <2e-16 ***
ORIGIN_SZNTSZ02 16.5792122 0.0233186 710.99 <2e-16 ***
ORIGIN_SZNTSZ03 18.9506415 0.0075957 2494.93 <2e-16 ***
ORIGIN_SZNTSZ05 15.8770261 0.0495825 320.21 <2e-16 ***
ORIGIN_SZNTSZ06 15.3997415 0.0557029 276.46 <2e-16 ***
ORIGIN_SZNVSZ01 20.2241694 0.0043487 4650.65 <2e-16 ***
ORIGIN_SZNVSZ02 19.1897826 0.0065383 2934.97 <2e-16 ***
ORIGIN_SZNVSZ03 18.8854268 0.0080459 2347.22 <2e-16 ***
ORIGIN_SZNVSZ04 18.8940191 0.0090985 2076.61 <2e-16 ***
ORIGIN_SZNVSZ05 17.6278585 0.0168107 1048.61 <2e-16 ***
ORIGIN_SZPGSZ01 19.4825220 0.0122960 1584.46 <2e-16 ***
ORIGIN_SZPGSZ02 19.4726761 0.0073116 2663.25 <2e-16 ***
ORIGIN_SZPGSZ03 20.5515713 0.0045631 4503.86 <2e-16 ***
ORIGIN_SZPGSZ04 21.0527131 0.0041500 5072.89 <2e-16 ***
ORIGIN_SZPGSZ05 20.1436604 0.0057267 3517.48 <2e-16 ***
ORIGIN_SZPLSZ01 19.1832002 0.0120006 1598.53 <2e-16 ***
ORIGIN_SZPLSZ02 18.8752206 0.0149740 1260.53 <2e-16 ***
ORIGIN_SZPLSZ03 18.1000818 0.0371769 486.86 <2e-16 ***
ORIGIN_SZPLSZ04 17.1730559 0.0370280 463.79 <2e-16 ***
ORIGIN_SZPLSZ05 17.9084439 0.0225031 795.82 <2e-16 ***
ORIGIN_SZPNSZ01 21.0804425 0.0044829 4702.41 <2e-16 ***
ORIGIN_SZPNSZ02 19.8822123 0.0111507 1783.05 <2e-16 ***
ORIGIN_SZPNSZ03 17.9293289 0.0193571 926.24 <2e-16 ***
ORIGIN_SZPNSZ04 17.1039594 0.0334954 510.64 <2e-16 ***
ORIGIN_SZPNSZ05 18.2543864 0.0275554 662.46 <2e-16 ***
ORIGIN_SZPRSZ01 19.8777935 0.0117586 1690.49 <2e-16 ***
ORIGIN_SZPRSZ02 21.0751780 0.0044832 4700.88 <2e-16 ***
ORIGIN_SZPRSZ03 20.6717019 0.0045577 4535.55 <2e-16 ***
ORIGIN_SZPRSZ04 19.6365125 0.0074923 2620.90 <2e-16 ***
ORIGIN_SZPRSZ05 21.3132151 0.0042119 5060.24 <2e-16 ***
ORIGIN_SZPRSZ06 18.9314574 0.0117278 1614.24 <2e-16 ***
ORIGIN_SZPRSZ07 17.2822918 0.0162430 1063.98 <2e-16 ***
ORIGIN_SZPRSZ08 19.9267642 0.0062298 3198.62 <2e-16 ***
ORIGIN_SZQTSZ01 19.7357175 0.0066359 2974.08 <2e-16 ***
ORIGIN_SZQTSZ02 19.2082141 0.0061402 3128.26 <2e-16 ***
ORIGIN_SZQTSZ03 19.7771883 0.0056220 3517.83 <2e-16 ***
ORIGIN_SZQTSZ04 18.7114421 0.0072842 2568.76 <2e-16 ***
ORIGIN_SZQTSZ05 19.3049324 0.0062401 3093.69 <2e-16 ***
ORIGIN_SZQTSZ06 19.2643228 0.0065590 2937.09 <2e-16 ***
ORIGIN_SZQTSZ07 18.5697347 0.0095373 1947.06 <2e-16 ***
ORIGIN_SZQTSZ08 19.6147001 0.0061330 3198.21 <2e-16 ***
ORIGIN_SZQTSZ09 19.2550793 0.0069947 2752.82 <2e-16 ***
ORIGIN_SZQTSZ10 19.5801866 0.0064513 3035.07 <2e-16 ***
ORIGIN_SZQTSZ11 17.7398366 0.0143648 1234.95 <2e-16 ***
ORIGIN_SZQTSZ12 17.2420354 0.0186736 923.34 <2e-16 ***
ORIGIN_SZQTSZ13 19.3857418 0.0078878 2457.69 <2e-16 ***
ORIGIN_SZQTSZ14 18.1300753 0.0122096 1484.90 <2e-16 ***
ORIGIN_SZQTSZ15 19.4222283 0.0120871 1606.86 <2e-16 ***
ORIGIN_SZRCSZ01 18.1549045 0.0125108 1451.13 <2e-16 ***
ORIGIN_SZRCSZ06 18.8836400 0.0082161 2298.38 <2e-16 ***
ORIGIN_SZRVSZ01 16.7864438 0.0323796 518.43 <2e-16 ***
ORIGIN_SZRVSZ02 16.4203244 0.0276836 593.14 <2e-16 ***
ORIGIN_SZRVSZ03 16.6453738 0.0244992 679.42 <2e-16 ***
ORIGIN_SZRVSZ04 15.9559213 0.0556344 286.80 <2e-16 ***
ORIGIN_SZRVSZ05 17.0476331 0.0164122 1038.71 <2e-16 ***
ORIGIN_SZSBSZ01 20.0417968 0.0062488 3207.29 <2e-16 ***
ORIGIN_SZSBSZ02 19.1869565 0.0081051 2367.26 <2e-16 ***
ORIGIN_SZSBSZ03 20.5769861 0.0045108 4561.70 <2e-16 ***
ORIGIN_SZSBSZ04 20.5154199 0.0050548 4058.57 <2e-16 ***
ORIGIN_SZSBSZ05 19.6250669 0.0065562 2993.35 <2e-16 ***
ORIGIN_SZSBSZ06 18.8419757 0.0171135 1101.00 <2e-16 ***
ORIGIN_SZSBSZ07 19.4897259 0.0124528 1565.09 <2e-16 ***
ORIGIN_SZSBSZ08 18.7027917 0.0140545 1330.73 <2e-16 ***
ORIGIN_SZSBSZ09 18.8893480 0.0088571 2132.67 <2e-16 ***
ORIGIN_SZSESZ02 20.8962192 0.0041665 5015.34 <2e-16 ***
ORIGIN_SZSESZ03 20.9452771 0.0039737 5270.94 <2e-16 ***
ORIGIN_SZSESZ04 20.6576142 0.0046364 4455.55 <2e-16 ***
ORIGIN_SZSESZ05 19.5170732 0.0058912 3312.92 <2e-16 ***
ORIGIN_SZSESZ06 20.7595824 0.0045747 4537.89 <2e-16 ***
ORIGIN_SZSESZ07 17.6888256 0.0195787 903.47 <2e-16 ***
ORIGIN_SZSGSZ01 19.1359250 0.0085781 2230.79 <2e-16 ***
ORIGIN_SZSGSZ02 18.5614369 0.0102037 1819.10 <2e-16 ***
ORIGIN_SZSGSZ03 19.9933176 0.0050434 3964.23 <2e-16 ***
ORIGIN_SZSGSZ04 20.2426871 0.0047211 4287.71 <2e-16 ***
ORIGIN_SZSGSZ05 18.0114965 0.0107743 1671.70 <2e-16 ***
ORIGIN_SZSGSZ06 20.2593194 0.0044538 4548.76 <2e-16 ***
ORIGIN_SZSGSZ07 19.0763664 0.0062968 3029.54 <2e-16 ***
ORIGIN_SZSKSZ01 19.9222451 0.0085136 2340.04 <2e-16 ***
ORIGIN_SZSKSZ02 20.8633383 0.0055248 3776.33 <2e-16 ***
ORIGIN_SZSKSZ03 19.6528148 0.0080534 2440.33 <2e-16 ***
ORIGIN_SZSKSZ04 18.0754470 0.0275771 655.45 <2e-16 ***
ORIGIN_SZSKSZ05 19.1192521 0.0155579 1228.91 <2e-16 ***
ORIGIN_SZSLSZ01 17.1501034 0.0329384 520.67 <2e-16 ***
ORIGIN_SZSLSZ04 19.5949774 0.0076753 2552.98 <2e-16 ***
ORIGIN_SZSRSZ01 16.9761403 0.0162020 1047.78 <2e-16 ***
ORIGIN_SZTHSZ01 17.9695687 0.0488559 367.81 <2e-16 ***
ORIGIN_SZTHSZ03 18.5427522 0.0223617 829.22 <2e-16 ***
ORIGIN_SZTHSZ04 17.4760374 0.0286247 610.52 <2e-16 ***
ORIGIN_SZTHSZ06 17.8401186 0.0183322 973.16 <2e-16 ***
ORIGIN_SZTMSZ01 20.3406361 0.0056607 3593.33 <2e-16 ***
ORIGIN_SZTMSZ02 22.0307026 0.0037386 5892.85 <2e-16 ***
ORIGIN_SZTMSZ03 21.3451920 0.0040606 5256.65 <2e-16 ***
ORIGIN_SZTMSZ04 20.6611593 0.0049896 4140.87 <2e-16 ***
ORIGIN_SZTMSZ05 19.3323133 0.0112868 1712.82 <2e-16 ***
ORIGIN_SZTNSZ01 17.9513571 0.0128266 1399.54 <2e-16 ***
ORIGIN_SZTNSZ02 18.0267387 0.0098372 1832.51 <2e-16 ***
ORIGIN_SZTNSZ03 17.7253700 0.0134668 1316.23 <2e-16 ***
ORIGIN_SZTNSZ04 19.4474075 0.0073760 2636.59 <2e-16 ***
ORIGIN_SZTPSZ01 19.1078631 0.0065635 2911.25 <2e-16 ***
ORIGIN_SZTPSZ02 20.2837634 0.0041411 4898.18 <2e-16 ***
ORIGIN_SZTPSZ03 19.1838238 0.0059552 3221.37 <2e-16 ***
ORIGIN_SZTPSZ04 19.1805388 0.0054778 3501.53 <2e-16 ***
ORIGIN_SZTPSZ05 19.3718076 0.0058610 3305.18 <2e-16 ***
ORIGIN_SZTPSZ06 19.6605723 0.0054968 3576.70 <2e-16 ***
ORIGIN_SZTPSZ07 19.4499807 0.0060491 3215.36 <2e-16 ***
ORIGIN_SZTPSZ08 18.7996538 0.0095757 1963.28 <2e-16 ***
ORIGIN_SZTPSZ09 19.0025110 0.0067068 2833.31 <2e-16 ***
ORIGIN_SZTPSZ10 18.8899657 0.0076094 2482.46 <2e-16 ***
ORIGIN_SZTPSZ11 19.6277780 0.0053983 3635.93 <2e-16 ***
ORIGIN_SZTPSZ12 19.1471104 0.0065742 2912.45 <2e-16 ***
ORIGIN_SZTSSZ01 17.4901113 0.0478954 365.17 <2e-16 ***
ORIGIN_SZTSSZ02 20.4997466 0.0081850 2504.55 <2e-16 ***
ORIGIN_SZTSSZ03 20.1076553 0.0084728 2373.19 <2e-16 ***
ORIGIN_SZTSSZ04 20.0646610 0.0089008 2254.26 <2e-16 ***
ORIGIN_SZTSSZ05 19.3962067 0.0151392 1281.19 <2e-16 ***
ORIGIN_SZTSSZ06 20.9235857 0.0178278 1173.65 <2e-16 ***
ORIGIN_SZWCSZ01 20.8411600 0.0086519 2408.86 <2e-16 ***
ORIGIN_SZWCSZ02 17.7355404 0.0328889 539.26 <2e-16 ***
ORIGIN_SZWCSZ03 14.9380886 0.1240699 120.40 <2e-16 ***
ORIGIN_SZWDSZ01 21.1969012 0.0037830 5603.23 <2e-16 ***
ORIGIN_SZWDSZ02 20.5930001 0.0044572 4620.13 <2e-16 ***
ORIGIN_SZWDSZ03 21.2521867 0.0041672 5099.85 <2e-16 ***
ORIGIN_SZWDSZ04 21.0702687 0.0048648 4331.13 <2e-16 ***
ORIGIN_SZWDSZ05 20.4008998 0.0051801 3938.35 <2e-16 ***
ORIGIN_SZWDSZ06 20.6669176 0.0049280 4193.78 <2e-16 ***
ORIGIN_SZWDSZ07 19.0500370 0.0082729 2302.71 <2e-16 ***
ORIGIN_SZWDSZ08 19.0816252 0.0080667 2365.49 <2e-16 ***
ORIGIN_SZWDSZ09 21.4182096 0.0040391 5302.73 <2e-16 ***
ORIGIN_SZYSSZ01 19.5355157 0.0057540 3395.14 <2e-16 ***
ORIGIN_SZYSSZ02 20.8737972 0.0048278 4323.64 <2e-16 ***
ORIGIN_SZYSSZ03 21.6614437 0.0040011 5413.81 <2e-16 ***
ORIGIN_SZYSSZ04 20.9305289 0.0043595 4801.10 <2e-16 ***
ORIGIN_SZYSSZ05 20.1727678 0.0058466 3450.34 <2e-16 ***
ORIGIN_SZYSSZ06 19.1481507 0.0116724 1640.47 <2e-16 ***
ORIGIN_SZYSSZ07 18.7919074 0.0141636 1326.78 <2e-16 ***
ORIGIN_SZYSSZ08 19.9733515 0.0061229 3262.07 <2e-16 ***
ORIGIN_SZYSSZ09 20.9366181 0.0040347 5189.15 <2e-16 ***
log(SCHOOL_COUNT) 0.4755516 0.0004701 1011.55 <2e-16 ***
log(RETAIL_COUNT) 0.1796905 0.0001856 968.12 <2e-16 ***
log(DIST) -1.6929522 0.0004093 -4136.01 <2e-16 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
(Dispersion parameter for poisson family taken to be 1)
Null deviance: 189463537 on 14471 degrees of freedom
Residual deviance: 15526121 on 14189 degrees of freedom
AIC: 15615824
Number of Fisher Scoring iterations: 6Check the log variables, 1st col should always be positive to attract people. Distance must be -ve bc it is inverse relationship. ie people less willing to travel with longer dist = distance decay p-values should be < 0.05 to accept as part of the conceptual models. if > 0.05 = not statistically significant.
8.1. Goodness-of-Fit
CalcRSquared <- function(observed, estimated) {
r <- cor(observed, estimated)
R2 <- r^2
R2
}
CalcRSquared(orcSIM_Poisson$data$TRIPS, orcSIM_Poisson$fitted.values)[1] 0.4362208Used residuals in the orcSIM_Poisson table to see which one is useful
8.2 Root mean square error
9. Doubly Constrained
Dont need to -1 because no destination attractiveness
10. Model Comparison
Find out outlier, see how model is affected if outlier is removed.